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ABSTRACT 

Using the sample of long Gamma-ray bursts (GRBs) detected by Swift-BAT 
before June 2007, we measure the cumulative distribution of the peak photon 
fluxes (log iV-log P) of the Swift bursts. Compared with the BATSE sample, 
we find that the two distributions are consistent after correcting the band pass 
difference, suggesting that the two instruments are sampling the same population 
of bursts. We also compare the log AMog P distributions for sub-samples of the 
Swift bursts, and find evidence for a deficit (99.75% confident) of dark bursts 
without optical counterparts at high peak flux levels, suggesting different red- 
shift or 7-ray luminosity distributions for these bursts. The consistency between 
the logiV-logP distributions for the optically detected bursts with and with- 
out redshift measurements indicates that the current sample of the Swift bursts 
with redshift measurements, although selected heterogeneously, represents a fair 
sample of the non-dark bursts. We calculate the luminosity functions of this 
sample in two redshift bins (z < 1 and z > 1), and find a broken power-law 
is needed to fit the low redshift bin, where dN/dL oc £- L27±006 for high lumi- 
nosities (Lp ea k > 5 x 10 48 erg s _1 ) and dN/dL oc £- 2 - 3±0 - 3 at for low luminosities, 
confirming the results of several studies for a population of low luminosity GRBs. 

Subject headings: gamma rays: bursts 

1. Introduction 

The cumulative distribution of source intensities (log iV-log S) is a useful tool for study- 
ing the source populations, especially when the redshifts of the sources are not measured. In 
the Gamma-ray burst (GRB) field, since the peak photon flux of the GRB is directly related 
to the detection threshold, the cumulative distribution of peak photon flux (log AMog P) is 
used in many studies. A number of these studies have made use of the BATSE sample of 
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more than 2,000 GRBs. This sample was first used to demonstrate the cosmological nature 
of the bursts (e.g., Fenimore et al. 1993; Pendleton et al. 1996), and then used to constrain 
the GRB populations (e.g., Kommers et al. 2000; Lin et al. 2004; Guetta et al. 2005; Dai & 
Zhang 2005), though many assumed that GRB rate follows the star formation rate in the 
analyses. Since the launch of Swift (Gehrels et al. 2004), the sample of the Swift-B AT bursts 
has reached to 237 in June 2007 (Sakamoto et al. 2008). With the sample size, it is now pos- 
sible to compare the Swift and BATSE samples to test whether the two instruments sample 
the same population of bursts. It is possible that they are different since the two instruments 
have different band passes and sensitivities. For example, Band (2006) showed that Swift 
is more sensitive to soft bursts. Moreover, there is still a large number of bursts without 
redshift measurements, e.g., dark bursts, and the logiV-logP distribution provides a means 
of studying their source population besides using pseudo-redshifts derived from spectral or 
timing properties (e.g., Norris 2002). 

The fraction of the Swift bursts with redshift measurements has increased significantly 
compared to the BATSE bursts. Although selected heterogeneously, it is tempting to mea- 
sure the luminosity function of the Swift bursts using this sample (e.g., Liang et al. 2007). 
Besides issues with redshift selection effects, the Swift trigger efficiency has not been well 
studied, which presents an additional difficulty. We show that by studying the logiV-logP 
distributions of the Swift and BATSE bursts, we can justify the usage of the heterogeneous 
redshift sample and set detection thresholds for measuring luminosity functions. We present 
the luminosity functions using the heterogeneous redshift sample, where we adopt a cosmol- 
ogy of H = 70 km s' 1 Mpc" 1 , fi m = 0.3, and fi A = 0.7. 



2. The Swift Burst Sample 

We use the Swift-B AT GRB catalog published in Sakomoto et al. (2008). The sample 
consists of 237 Swift bursts detected before 2007 June 16. All the bursts are triggered by 
the BAT instrument on board Swift. The Swift-BAT catalog contains a number of basic 
properties of the bursts such as the burst duration, spectral index, and peak photon flux 
in several bands. There are 229 bursts with peak photon flux estimates in the 15-150 keV 
band, with a minimum value of 0.23 photon cm~ 2 s -1 . Of these 229 bursts, 210 bursts can 
be identified as long bursts and 15 as short bursts, where we divide the sample at T 90 = 2 s. 
In this paper, we focus on these 210 long bursts. To compare with the BATSE burst sample, 
we use the long BATSE burst sample from Kommers et al. (2000), which consists of 2176 
long GRBs from both online and off-line searches. 
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Fig. 1. — The logAMogP distributions of the Swift bursts (thick solid line) and BATSE 
bursts (thin solid line). We use the Swift sample from Sakamoto et al. (2008) and the 
BATSE sample from Kommers et al. (2000). The peak flux values of the BATSE bursts 
are corrected to the BAT band pass of 15-150 keV using an average correction factor of 
A50-300 = 0.4 * Ai5_i5 . The K-S test shows that the null probability that the two samples 
are drawn from the same distribution is 69%, which indicates that the BAT and BATSE 
instruments sample the same population of GRBs. 
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3. The log A^— log P Distribution 

We plot the log AMog P distributions of the Swift and BATSE bursts in Figure [TJ We 
use the > 50% coded field of view of ^ Swift = s ^ r f° r the BAT and a flight time of 
T Swift = ^ r un til 2007 June 16 to calculate the burst rate for the Swift sample. For the 
BATSE sample, we adopt the values from Kommers et al. (2000) with Tbatse — 1-33 x 10 8 s 
and a mean field of view of flsATSE = 0.67 * 4n str. We also correct for the band pass 
differences between the BAT (15-150 keV) and BATSE (50-300 keV) instruments, where 
we use the spectral fits provided by Sakamoto et al. (2008) for the Swift sample. Using the 
simple power-law fits to the BAT spectra, which fit well for most Swift bursts (Sakamoto 
et al. 2008), we find a mean relation, A50-300 = 0.4 * Ni 5 _i 50 , which we use to correct the 
photon flux for the BATSE bursts. Figure [T] shows that the two distributions are quite 
consistent, and only small discrepancies with ~ la significance exist at high photon flux 
levels at f pea k > 8 photon cm -2 s _1 . We perform a Kolmogorov-Smirnov (K-S) test for 
the two samples, and find that the null probability that the two samples are drawn from 
the same distribution is 69% (76% for the sub-samples with f vea k > 8 photon cm -2 s _1 ). 
This result indicates that the BAT and BATSE instruments sample the same population of 
GRBs, even though they have different band passes. We also convert the BAT photon flux 
to the BATSE band individually for each burst using the simple power-law fits of Sakamoto 
et al. (2008) to reduce the uncertainty introduced by scatter of the spectral indices around 
their mean, and compare the two distributions in the BATSE band. We again find that 
the two distributions are not significantly different. For some bursts, the spectra are better 
fit by a cut-off power- law (Sakamoto et al. 2008), which introduces additional uncertainties 
in the flux conversion; however they are generally negligible compared to the uncertainty 
introduced by scatter of the photon indices. We note that the detection limits of the two 
samples also match each other after the band pass correction. 



3.1. Comparisons between Swift Sub-Samples 

Next, we compare the log AMog P distributions for sub-samples of the Swift bursts. 
First, we compare the first half and the second half of the Swift bursts separated by time, 
finding that the two samples are consistent, suggesting that there is no significant degrada- 
tion in the BAT detector. Next, we compare the distributions for bursts with or without 
redshift measurements or optical detections (Figure [2]). There are 68 bursts with redshift 
measurements and optical afterglow detections, 48 with optical afterglow detections but with- 
out redshift measurements, and 94 without optical afterglow detections. The last category 
of bursts with no optical afterglow detections is the dark bursts, though another definition 
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Fig. 2. — Comparisons between the log TV-log P distributions for the Swift sub-samples. The 
thick solid line is for the Swift bursts with optical detections and redshift measurements, and 
the thick dashed line is for bursts with optical detections but without redshift measurements. 
The thin dotted line is for the dark bursts without optical detections, and the thin dash- 
dotted line is for bursts with UVOT observations and with Galactic extinction Ay t Gai < 
1.5 mag, but without optical detections. Although the K-S test results show that all four 
sub-samples are consistent with K-S probabilities of ~ 35%, there is an indication that the 
dark bursts are under-populated at high peak flux levels. The K-S test probability drops 
to 4% when we compare the high peak flux bursts (f pea k > 5 photon cm -2 s _1 ) between 
the distributions for dark bursts and optically detected bursts. At even higher peak flux 
levels (fpeak > 10 photon cm~ 2 s _1 ), the difference is more pronounced with a Poisson null 
probability of 0.25%. 
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uses an optical to X-ray threshold ratio to define this sample (Jakobsson et al. 2004). In this 
paper, we use the simple definition of no optical detection for dark bursts. In most cases, 
we find no significant difference between the sub-samples based on the K-S test probabilities 
(45% between samples with and without redshift measurements, 34% between samples with 
redshift measurements and dark bursts, and 24% between samples with redshift measure- 
ments and without redshift measurements but with optical detections). 

However, Figure [2] shows that there is an indication for a deficit of dark bursts at high 
peak flux levels compared to the bursts with optical afterglow detections. We test it by 
including a flux filter and select bursts with f pea k > 5 photon cm -2 s" 1 , and find that the K- 
S probability that the dark bursts and bursts with optical detections are drawn from the same 
population is 4%. It seems that the deficit is more pronounced in the flux level at f pea k > 
10 photon cm -2 s _1 . Since the number of bursts detected in this regime is too small and 
that the K-S test is no longer applicable, we use a simple Poisson argument instead. The null 
model is that the dark bursts follow the same distribution as the optically detected bursts, 
which predicts that the number of the dark bursts detected at f pea k > 10 photon cm -2 s _1 
is 8.2±0.8. This model prediction and the associated error-bar are obtained by aligning the 
normalizations of the two distributions at low flux levels. Since we only detect one dark 
burst with f pea k > 10 photon cm -2 s _1 , the Poisson probability of detecting no more than 
one burst is 0.0025 for an expected value of 8.2. If we use the lower end of the model 
prediction (7.4), the corresponding probability is 0.0051. We argue that there is evidence 
that the dark bursts do not follow the logiV-logP distribution of the optically detected 
bursts. It can be either a deficit of dark bursts at high peak flux levels or an over abundance 
of dark bursts at low flux levels. Future analysis including the whole Swift sample is needed 
to confirm this results. 

Since the ground-based optical follow-up observations of GRBs can be potentially biased 
by several reasons, such as the scheduling issues or unfavorable burst sky positions, we also 
test by limiting the dark burst sample using those bursts observed by UVOT but without 
UVOT or other optical detections (UVOT dark bursts). This provides a more uniform 
selection because UVOT routinely observes the GRB fields after the burst triggers. In 
addition, this excludes some of the bursts with unfavorable sky positions to follow-up in 
the optical bands, such as those located close to the Sun. We also limit the sample by 
excluding bursts with large Galactic extinctions, Ay,Gai > 1-5 mag, where the dark burst 
fraction is significantly larger compared to that for the total population. However, since the 
UVOT flux limit is generally shallower than those of the ground-based observations using 
mid-sized telescopes, we still classify the dark bursts as those without optical detections. 
We find the logiV-logP distribution of the UVOT dark burst sample follows that of the 
general dark burst sample (or no optical counterparts sample, Figure [2]). This is because 
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only a small fraction of bursts have unfavorable burst positions that cannot be observed by 
the ground-based telescopes. In addition, a number of optical telescopes are used to observe 
GRB afterglows world wide, which collectively reduces the scheduling problems such as the 
whether conditions in individual optical sites. Because of the shallowness of the UVOT flux 
limit, it is still possible that both of the dark burst (no optical counterpart) samples contain 
a fraction of normal, non-dark bursts, because they are not observed promptly by the mid 
or large size ground-based telescopes. However, this population should not contribute to the 
difference between the log iV-log P distributions between dark and normal bursts. 



4. Luminosity Function 

The bursts with redshift measurements account for 59% of the bursts with optical af- 
terglow detections. Figure [2] shows that the log iV-log P distributions between the optically 
detected bursts with and without redshift measurements are quite similar. This suggests that 
the current redshift sample, although obtained heterogeneously, represents a fair sample of 
the bursts with optical afterglow detections. We compute the luminosity function (LF) of this 
sample using the l/V max method (Schmidt 1968; Avni & Bahcall 1980), where we include an 
additional factor of 1/(1 +z) in the differential volume element to account for the cosmo logical 
time dilation effect. We calculate the k-corrections using the power-law fit to the BAT spec- 
tra from Sakomoto et al. (2008), and use a detection limit of f pea k = 0.25 photon cm -2 s _1 , 
the smallest value in the redshift sample (GRB 060218), to calculate the V max values. Since 
the log iV-log P distribution of the Swift sample matches that of the BATSE sample when 
approaching the detection limits (Figured]), we use the trigger efficiency analysis of Kom- 
mers et al. (2000) to estimate the efficiency of the Swift triggers, and then assign weights to 
the Swift bursts with a maximum weight of 2. Using a maximum weight is a conservative 
approach to avoid huge corrections at low peak flux levels, where the uncertainties of the 
efficiency analysis could be large. 

We show the total LF and those in redshift bins of z < 1 and z > 1 in Figure [31 where 
we choose a bin size of 0.5 dex. We find that the LF for the lower redshift bin is consistent 
with the total LF, and they cannot be fit well by a single power-law model with x 2 /dof = 
4.6 and 2.8, respectively, for the total and z < 1 LFs. We add another power-law component 
at the low luminosity end, and fit the z < 1 LF. We find dN/dL oc £™ 2 - 3±a3 at the low 
luminosity end and dN/dL oc /,- 1 - 27±0 06 a t the high luminosity end {L vea k > 5 x 10 48 erg s _1 ) 
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with x 2 /dof =1.3. The complete LF for the z < 1 bin is in Equation [TJ 

r \ -2.3±0.3 / L \ _L27±a06 ~ 

5 x 10 48 erg s^ 1 J + \5 x 10 48 erg s" 1 / 

(1) 

Comparing the LFs in the two redshift bins, the two low luminosity data points for the z > 1 
LF at L pea /c = 10 50_51 erg s _1 are below the z < 1 LF. Since the star formation rate drops 
significantly below z = 1 (e.g., Hopkins & Beacom 2006), if the GRB rate follows the star 
formation rate, we expect that the LF for the z > 1 bin should have a higher normalization 
than that for the z < 1 LF. It is possible that the two low luminosity data points of the 
z > 1 LF are incomplete because they are close to the Swift flux limit (e.g., Kistler et al. 
2008), even though we have partially modeled the Swift trigger incompleteness by using 
the BATSE trigger efficiency model (Kommers et al. 2000). If we neglect these two low 
luminosity data points, the remaining two data points of the z > 1 LF are not lower than 
the z < 1 LF. Considering the uncertainties, the results do not present a significant challenge 
to the hypothesis that the GRB rate follows the star formation rate. We note that the LFs 
do not include the contribution from the dark bursts. We need to multiply by a factor of 
1.8 to include them; however, this may not be accurate since the dark bursts may have a 
different redshift or 7-ray luminosity distribution (§3.1). In addition, if the majority of the 
dark bursts are at high redshifts (§5), the normalization of the z > 1 LF will be significantly 
increased compared to the z < 1 LF, which can raise the GRB rate to be consistent with 
the star formation rate, or even evolving faster as suggested by some recent studies (e.g., 
Salvaterra & Chinarini 2007; Kistler et al. 2008; Salvaterra et al. 2009). 



^ = (1.7±0.2)l4Gpc- 3 Yr- 1 



5. Discussion 

We compare the log AMog P distributions for the Swift bursts and BATSE bursts, and 
find that they are consistent after correcting the band pass differences as suggested by the 
K-S test results. This shows that the two instruments sample the same population of bursts, 
although they have different band passes. We can also compare the normalizations of the 
distributions directly, since the field of view and the observing time of the two samples are 
measured, and we find that they are consistent. Indirect comparison between the Swift 
and BATSE samples has been performed in previous studies (Salvaterra & Chincarini 2007; 
Virgili et al. 2009), where the two samples are found to be consistent with same theoretical 
models. Given that Swift is more sensitive than BATSE for the soft bursts (Band 2006), 
it is puzzling that the two distributions are consistent. Currently, we are comparing the 
on-board triggered Swift sample with the BATSE off-line search sample. It is possible that 
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Fig. 3. — The luminosity function of the Swift bursts with redshift measurements. The data 
for the z < 1 and z > 1 bins are slightly shifted for clarity. The dashed dotted line is a 
broken power-law fit to the z < 1 LF. 
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the difference will become significant when comparing both the off-line search samples. 

We also compare the log TV-log P distributions for the sub-samples of the Swift bursts 
with or without optical afterglow detections or redshift measurements. We find them to be 
broadly consistent. However, at the high peak flux regime (f pea k > 5 photon cm~ 2 s -1 ), we 
find an indication that there is a deficit (96% confident) of "dark bursts" (bursts without 
optical counterparts). This deficit is more evident at f pea k > 10 photon cm~ 2 s -1 , where we 
find the effect is significant at the 99.75% confidence level. 

A comparison between the log AMog P distributions for dark bursts and optically bright 
bursts can place important constraints on the origin of the dark bursts. There are several 
models proposed for the dark bursts. They can be bursts with intrinsic normal optical-to-7- 
ray ratios but with large intrinsic or foreground optical extinction (e.g., Taylor et al. 1998; 
Djorgovski et al. 2001; Fynbo et al. 2001), have intrinsic faint optical afterglows, i.e., low 
optical-to-7-ray ratios, (e.g., Groot et al. 1998a; Frail et al. 1999), are normal bursts but exist 
in extremely high redshifts where the Lyman break lands in the optical (e.g., Groot et al. 
1998b; Fynbo et al. 2001), or a combination of several origins mentioned above. A different 
7-ray luminosity distribution or redshift distribution is needed to interpret the difference 
in the log AMog P distributions. Since the 7-ray/hard X-ray flux is not sensitive to the 
ISM absorption, the hypotheses with optical extinction or intrinsic optical faintness predict 
that the log AMog P distributions should have a similar shape between dark and optically 
detected bursts. For the high redshift scenario, it is possible for the dark bursts to have a 
different shape in the log AMog P distribution from the normal, non-dark bursts. Therefore, 
the deficit of the dark bursts at high flux levels suggests that a significant fraction of dark 
burst is from extreme high redshifts. In addition, it is possible that the bulk of the dark 
bursts are selected due to their low optical luminosity, and the deficit of dark bursts at high 
peak photon fluxes merely reflects the detection threshold of optical observations under the 
assumption of a constant optical-to-7-ray ratio. 

We measure the luminosity function of the Swift bursts using the sample with redshift 
measurements. We find the luminosity function can be fit by a broken power-law with 
dN/dL oc L- 2 - 3±0 - 3 at the low luminosity end, and dN/dL oc L _1 - 27±0 - 06 at the high luminosity 
end. We compare this result with previous measurements (Schmidt 2001; Norris 2002; Stern 
et al. 2002; Firmani et al. 2004; Guetta et al. 2005; Liang et al. 2007; Salvaterra & Chincarini 
2007) and find a best match with the result from Liang et al. (2007). The requirement for 
an additional component at the low luminosity end confirms the existence of the population 
of low luminosity GRBs claimed by several studies (e.g., Cobb et al. 2006; Piran et al. 2006; 
Soderberg et al. 2006; Guetta & Delia Valle 2007; Liang et al. 2007; Virgili et al. 2009). 
Although we model it as a power-law, it could be the tail of a Gaussian component. At the 
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high luminosity range, the measured slope dN/dL oc £- L27 ± - 06 i s close to the prediction 
of the "quasi-universal Gaussian jet" (dN/dL oc L -1 , Zhang et al. 2004). At the very high 
luminosity end (L peak > 10 51 erg s -1 ), the "quasi-universal Gaussian jet" predicts dN/dL oc 
L~ 2 (Lloyd- Ronning et al. 2004; Dai & Zhang 2004), which cannot be tested with the current 
data. 

Unlike many of the previous studies, we do not make any assumption on the GRB rate, 
since there are recent studies suggesting that the GRB rate does not follow the star formation 
rate (e.g., Stanek et al. 2006). Instead, we measure the average GRB luminosity function in 
a large redshift bin. We compare the LFs in two luminosity bins. The z > 1 LF shows a 
drop at two low luminosity data points. This result can be affected by the uncertainties in 
the Swift trigger efficiency close to the Swift detection limit. If we neglect these two data 
points, the LFs for the two redshift bins are consistent. However, the large measurement 
uncertainties in the z < 1 LF make it difficult to test whether the GRB rate follows the star 
formation rate. 

We thank B. Zhang, C. S. Kochanek, R. Salvaterra, and the anonymous referee for 
helpful discussion. 

REFERENCES 

Avni, Y., & Bahcall, J. N. 1980, ApJ, 235, 694 
Band, D. L. 2006, ApJ, 644, 378 

Cobb, B. E., Bailyn, C. D., et al. 2006, ApJ, 645, L113 
Dai, X., & Zhang, B. 2005, ApJ, 621, 875 

Djorgovski, S. G., Frail, D. A., Kulkarni, S. R., et al. 2001, ApJ, 562, 654 

Fenimore, E. E., et al. 1993, Nature, 366, 40 

Firmani, C, Avila-Reese, V., et al. 2004, ApJ, 611, 1033 

Frail, D. A., et al. 1999, ApJ, 525, L81 

Fynbo, J. U., et al. 2001, A&A, 369, 373 

Gehrels, N., et al. 2004, ApJ, 611, 1005 



- 12 - 



Groot, P. J., et al. 1998a, ApJ, 502, L123 
Groot, P. J., et al. 1998b, ApJ, 493, L27 

Guetta, D., Granot, J., & Begelman, M. C. 2005, ApJ, 622, 482 

Guetta, D., & Delia Valle, M. 2007, ApJ, 657, L73 

Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 

Jakobsson, P., Hjorth, J., Fynbo, J. P. U., et al. 2004, ApJ, 617, L21 

Kommers, J. M., et al. 2000, ApJ, 533, 696 

Kistler, M. D., Yiiksel, H., et al. 2008, ApJ, 673, LI 19 

Liang, E., Zhang, B., Virgili, F., & Dai, Z. G. 2007, ApJ, 662, 1111 

Lin, J. R., Zhang, S. N., & Li, T. P. 2004, ApJ, 605, 819 

Lloyd-Ronning, N. M., Dai, X., & Zhang, B. 2004, ApJ, 601, 371 

Norris, J. P. 2002, ApJ, 579, 386 

Pendleton, G. N., et al. 1996, ApJ, 464, 606 

Pian, E., et al. 2006, Nature, 442, 1011 

Sakamoto, T., et al. 2008, ApJS, 175, 179 

Salvaterra, R., & Chincarini, G. 2007, ApJ, 656, L49 

Salvaterra, R., Guidorzi, C, et al. 2009, MNRAS, in press. larXiv:0805.4104l 

Schmidt, M. 1968, ApJ, 151, 393 

Schmidt, M. 2001, ApJ, 552, 36 

Soderberg, A. M., et al. 2006, Nature, 442, 1014 

Stanek, K. Z., et al. 2006, Acta Astronomica, 56, 333 

Stern, B. E., Tikhomirova, Y., & Svensson, R. 2002, ApJ, 573, 75 

Taylor, G. B., Frail, D. A., Kulkarni, S. R., et al. 1998, ApJ, 502, L115 

Virgili, F. J., Liang, E.-W., & Zhang, B. 2009, MNRAS, 392, 91 



-13- 

Zhang, B., Dai, X., et al. 2004, ApJ, 601, LI 19 



This preprint was prepared with the AAS IATgX macros v5.2. 



